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This  work  presents  time-averaged  ion  velocity  and  energy  distributions  in  the  axial  di¬ 
rection  within  the  acceleration  channel  and  axial  and  radial  directions  in  the  plume  of  a 
Princeton  University  low  power  cylindrical  Hall  thruster  operating  on  xenon.  Xenon  ion 
velocities  for  the  thruster  are  derived  from  laser-induced  fluorescence  measurements  of  the 
5d[4]7/2  —  6p[3]5/2  xenon  ion  excited  state  transition  at  A  =  843.72  nm.  Two  operating 
conditions  are  considered  with  variations  to  the  magnetic  field  strength,  in  an  effort  to 
capture  the  effect  on  ion  velocity  distributions.  The  lower  magnetic  field  condition  is  also 
considered  with  a  higher  vacuum  chamber  background  pressure.  Under  nominal  condi¬ 
tions,  xenon  ions  are  accelerated  to  a  most  probable  energy  of  25  eV  within  the  thruster 
with  an  additional  175  eV  gain  in  the  thruster  plume.  At  a  position  40  mm  into  the  plume, 
this  constitutes  an  energy  of  200  eV,  with  the  wings  of  the  energy  distribution  extending 
between  177  and  228  eV  at  an  applied  potential  of  300  V.  Decreasing  the  magnetic  field 
strength  appears  to  improve  axial  acceleration.  Increasing  the  background  chamber  pres¬ 
sure  pushes  the  region  of  maximum  acceleration  upstream  towards  the  thruster  exit  plane. 


I.  Introduction 

This  study  characterizes  time-averaged  xenon  ion  velocity  and  energy  distributions  both  inside  the  accel¬ 
eration  channel  and  in  the  near  field  plume  of  a  low  power  cylindrical  Hall  thruster  (CHT)  using  laser-induced 
fluorescence  (LIF).  It  is  a  continuation  of  a  previous  study  in  which  the  most  probable  axial,  radial  and  az¬ 
imuthal  velocities  within  the  acceleration  channel  and  in  the  plume  of  a  low  power  cylindrical  Hall  thruster 
were  reported.1  LIF  measurements  are  also  used  to  assess  the  effect  of  varying  magnetic  field  strength  on 
the  evolution  of  the  ion  velocity  and  kinetic  energy  profiles. 

The  electrostatic  thruster  used  in  this  study  is  the  Princeton  University  CHT-30.  Previous  efforts  exam¬ 
ining  low  power  annular  Hall  effect  thrusters  have  shown  that  erosion  of  the  insulator  protecting  the  central 
magnetic  pole  limits  thruster  lifetime  significantly,  even  more  so  than  erosion  of  the  outer  insulator.  This 
effect  is  amplified  as  the  power  and  size  of  the  Hall  thruster  is  scaled  down.  By  reducing  the  length  of  the 
central  magnetic  pole  and  creating  a  region  of  the  channel  with  a  low  surface-to- volume  area  (a  cylindrical 
region),  the  CHT  as  developed  by  Princeton  University  reduces  the  impact  of  insulator  erosion  significantly.2 

The  CHT-30  has  a  cylindrical  acceleration  channel  with  a  diameter  of  30  mm  and  depth  of  13  mm.  This 
is  similar  to  most  modern  Hall  thrusters,  which  have  acceleration  channels  with  a  maximum  depth  of  10-20 
mm.  With  this  geometry,  it  is  possible  to  align  collection  optics  to  the  probe  beam  such  that  limited  internal 

1  of  20 


American  Institute  of  Aeronautics  and  Astronautics 
Approved  for  public  release;  distribution  is  unlimited. 


optical  access  is  possible  without  thruster  modification.  In  this  work,  the  collection  lens  is  placed  60°  off 
the  plume  axis.3  With  the  CHT’s  plume  divergence  half  angle  of  approximately  60°  to  80°, 4  this  places 
the  collection  optics  away  from  much  of  the  plume  ion  flux.  In  this  way,  it  is  possible  to  probe  internal  ion 
acceleration  of  the  CHT-30  with  minimal  intrusion  into  the  plume. 

Two  thruster  operating  conditions  are  examined  in  this  study.  These  include  an  nominal  condition  for 
which  there  is  unpublished  performance  data,  and  an  reduced  magnetic  field  condition.  The  goals  of  these 
measurements  are  to  characterize  the  ion  velocity  for  a  known  case,  and  examine  a  low  plasma  oscillation 
case  where  the  signal-to-noise  ratio  (SNR)  of  the  diagnostic  is  maximized. 


II.  Experimental  Apparatus 


A.  Xenon  Ion  Spectroscopy 

For  the  results  reported  here,  the  5d[4]7/2  —  6p[3]5/2  electronic  transition  of  Xe  II  at  834.72  nm  is  probed.  The 
isotopic  and  nuclear-spin  effects  contributing  to  the  hyperfine  structure  of  the  5d[4]7/2  —  6p[3]5/2  xenon  ion 
transition  produce  a  total  of  19  isotopic  and  spin  split  components.  The  hyperfine  splitting  constants  which 
characterize  the  variations  in  state  energies  are  only  known  for  a  limited  set  of  energy  levels.  Unfortunately, 
this  xenon  ion  transition  only  has  confirmed  data  on  the  nuclear  spin  splitting  constants  of  the  6p[3]5/2  upper 
state.5-8 

Manzella  first  used  the  5d[4]7/2  —  6p[ 3] 5 /2  xenon  ion  transition  at  834.72  nm  to  make  velocity  measure¬ 
ments  in  a  Hall  thruster  plume.9  A  convenient  feature  of  this  transition  is  the  presence  of  a  relatively  strong 
line  originating  from  the  same  upper  state  (6s[2]3/2  —  6_p[3] 5 /2  transition  at  541.9  nm,10  which  allows  for 
non-resonant  fluorescence  collection) .  Ion  velocity  is  simply  determined  by  measurement  of  the  Doppler  shift 
of  the  absorbing  ions.11 


B.  Test  Facility 

The  LIF  measurements  were  performed  in  Chamber  6  of  the  Air  Force  Research  Laboratory  (AFRL)  Electric 
Propulsion  Laboratory  at  Edwards  AFB,  CA.  This  experimental  apparatus  has  been  described  extensively 
in  previous  work.1,3,12  Chamber  6  is  a  non-magnetic  stainless  steel  chamber  with  a  1.8  m  diameter  and 
3  m  length.  It  has  a  measured  pumping  speed  of  32,000  L/s  on  xenon.  Pumping  is  provided  by  four  single 
stage  cryogenic  panels  (single  stage  cold  heads  at  25  K)  and  one  50  cm  two  stage  cryogenic  pump  (12  K). 
Chamber  pressure  during  thruster  operation  is  approximately  2.8  x  10-6  torr,  corrected  for  xenon. 

Figure  1  shows  a  top  view  of  the  laser  optical  train,  collection  optics,  and  one  leg  of  the  external  probe 
optics.  The  laser  is  a  tunable  diode  laser,  capable  of  tuning  approximately  ±50  GHz  about  a  center  wave¬ 
length  of  834.72  nm.  The  laser  beam  passes  through  several  beam  pick-offs  followed  by  a  50-50  cube  beam 
splitter  where  it  is  split  into  two  beams  of  equal  power.  The  first  beam,  the  axial  probe  beam  shown  in 
Figures  1  and  2  is  focused  by  a  lens  and  enters  the  vacuum  chamber  through  a  window.  A  second  probe 
beam,  shown  in  Figure  2  only,  is  directed  from  the  optical  bench  via  a  periscope  apparatus  so  that  it  enters 
the  chamber  from  above  the  thruster  and  probes  the  velocity  perpendicular  to  the  first  probe  beam.  Each 
probe  beam  is  chopped  at  a  unique  frequency  by  choppers  Ch2  (2kHz)  and  Ch3  (2.8  kHz)  for  phase  sensitive 
detection  of  the  fluorescence  signals. 

The  two  wedge  beam  pick-offs  (BS)  shown  in  Figure  1  provide  portions  of  the  beam  for  diagnostic 
purposes.  These  diagnostics  include  a  commercial  wavelength  meter  and  a  300  MHz  free  spectral  range 
Fabry-Perot  etalon  (F-P)  to  monitor  wavelength  during  a  laser  scan,  and  a  low  pressure  xenon  hollow 
cathode  discharge  lamp  to  provide  a  stationary  absorption  reference  for  the  determination  of  the  Doppler 
shift.  Unfortunately,  there  is  no  detectable  population  of  the  ionic  xenon  5d[4]7/2  state.  However,  there  is  a 
nearby  (estimated  to  be  18.1  GHz  distant)  neutral  xenon  65' [1/2]  1  —  6p'[3/2]2  transition  at  834.68  nm. 13,14 

The  fluorescence  collection  optics  are  also  shown  in  Figure  1.  The  fluorescence  is  collected  by  a  75  mm 
diameter,  300  mm  focal  length  lens  within  the  chamber  and  oriented  60°  from  the  probe  beam  axis.  The 
collimated  fluorescence  signal  is  directed  through  a  window  in  the  chamber  side  wall  to  a  similar  lens  that 
focuses  the  collected  fluorescence  onto  the  entrance  slit  of  the  125  mm  focal  length  monochrometer  with  a 
photomultiplier  tube  (PMT).  Due  to  the  1:1  magnification  of  the  collection  optics,  the  spatial  resolution  of 
the  measurements  is  determined  by  the  geometry  of  the  entrance  slit  0.7  mm  width  and  1.5  mm  height  as 
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Figure  1.  Top  view  diagram  of  the  laser  optical  train  and  collection  optics.  Note  that  the  radial  probe  beam 
periscope  and  focusing  optics  are  not  shown. 


Figure  2.  Side  view  diagram  of  thruster  within  AFRL  chamber  6.  Also  shown  are  the  translation  stages  and 
the  laser  probe  beams.  Note  that  the  fluorescence  collection  and  external  optics  are  not  shown. 


3  of  20 


American  Institute  of  Aeronautics  and  Astronautics 
Approved  for  public  release;  distribution  is  unlimited. 


well  as  the  sub-mm  diameter  of  the  probe  beam.  This  apparatus  allows  for  limited  probing  of  the  interior 
acceleration  channel  of  Hall  thrusters  with  relatively  shallow  acceleration  channels.  Measurements  suggest 
that  this  combination  of  apparatus  and  laser  power  are  well  within  the  linear  fluorescence  regime. 

The  uncertainty  of  the  velocity  measurements  is  estimated  to  be  within  ±500  m/s.  The  repeatability 
of  the  peak  locations  appear  to  be  a  fraction  of  the  quoted  uncertainty.  However,  the  fluorescence  line 
shapes  are  often  significantly  broadened,  presumably  due  to  wide  velocity  distributions  caused  by  plasma 
fluctuations.  The  quoted  uncertainty  should  therefore  be  viewed  as  the  uncertainty  in  the  determination  of 
the  peak  of  the  fluorescence  line  shape. 

C.  Cylindrical  Hall  Thruster 

The  thruster  used  in  this  study  is  the  CHT-30  cylindrical  Hall  thruster  developed  by  Princeton  University. 
The  CHT-30  has  a  30  mm  diameter  channel  that  is  composed  of  a  short  annular  region  near  the  anode 
followed  by  a  13  mm  long  cylindrical  region  as  shown  in  Figure  3. 15  This  figure  also  indicates  the  magnetic 
field  topology  seen  during  the  operation  of  the  thruster.  CHTs  have  been  of  interest  in  the  electric  propulsion 
community  in  recent  years  due  to  their  low  power  operation  and  potential  for  long  life.  The  low  surface-to- 
volume  ratio  of  their  acceleration  channel  and  cusped  magnetic  field  configuration  result  in  low  erosion  rates 
of  their  ceramic  wall  materials,  which  is  especially  desirable  when  scaling  down  in  power.  CHTs  have  also 
demonstrated  performance  levels  approaching  those  of  state-of-the-art  annular  Hall  thrusters.2 


Figure  3.  CHT  Magnetic  field  distribution  for  cusped  field  configuration.15 


Although  the  CHT  and  several  other  magnetic  cusp  electrostatic  thrusters  have  elements  that  superficially 
resemble  Hall  effect  thrusters,  this  thruster  type  has  several  features  that  make  it  considerably  different.16 
First,  the  high  magnetic  field  strengths  within  cusp  thrusters  suggest  that  the  ions  may  be  magnetized.  Ion 
magnetization  is  understood  not  to  exist  in  annular  Hall  thrusters.  Furthermore,  the  geometry  and  magnetic 
field  structure  of  cusp  thrusters  are  notably  different,  with  the  cusped  magnetic  field  shape  seeing  no  analog 
in  traditional  Hall  thrusters.  In  this  early  stage  of  research,  these  differences  appear  to  indicate  that  the 
electron  transport  mechanism  within  cusp  thrusters  also  differs  significantly.  The  azimuthal  E  x  B  drift 
seen  in  conventional  Hall  thrusters  is  either  lacking  or  occurs  very  near  the  anode  of  the  cusp  thruster.  This 
suggests  that  the  E  x  B  shear  that  appears  critical  in  annular  Hall  thruster  electron  transport  does  not  exist, 
or  assumes  a  very  different  form  in  cusp  thrusters.1,2 

D.  Operating  Conditions 

The  operating  conditions  considered  in  this  study  were  a  nominal  operating  condition  provided  in  Table  1, 
as  well  as  an  reduced  magnetic  field  condition  provided  in  Table  2.  The  reduced  magnetic  field  operating 
condition  was  chosen  due  to  the  quiescent  operating  mode  of  the  thruster  and  the  improved  fluorescence 
signal  quality  throughout  the  majority  of  the  plume.  It  also  gave  insight  into  the  influence  of  the  magnetic 
field  strength  on  the  ion  velocity  profiles.  These  conditions  were  run  with  a  chamber  background  pressure  of 
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2.8  x  10  6  torr,  corrected  for  xenon.  The  reduced  magnetic  field  condition  was  also  run  at  a  higher  chamber 
background  pressure  of  3.5  x  10-5  torr  to  see  the  effect  of  varying  the  vacuum  chamber  conditions. 

Table  1.  Nominal  CHT-30  operating  conditions. 


Anode  Flow 

390  /ig/s  Xe 

Cathode  Flow 

98  fi g/s  Xe 

Anode  Potential 

300  V 

Anode  Current 

0.53  A 

Keeper  Current 

OFF 

Heater  Current 

5.0  A 

Magnet  1  Current 

3.00  A 

Magnet  2  Current 

2.00  A 

Table  2.  Reduced  Magnetic  Field  CHT-30  operating  conditions. 


Anode  Flow 

390  /ig/s  Xe 

Cathode  Flow 

150  /ig/s  Xe 

Anode  Potential 

300  V 

Anode  Current 

0.61  A 

Keeper  Current 

0.30  A 

Heater  Current 

5.0  A 

Magnet  1  Current 

2.00  A 

Magnet  2  Current 

2.00  A 

Figure  4  depicts  the  Cartesian  coordinate  system  used  in  these  measurements.  Figure  5  shows  a  side  view 
of  the  thruster  channel  and  plume  including  the  positions  of  the  LIF  data  points  examined  in  this  paper. 
The  thruster  is  mounted  on  a  three  axes  orthogonal  computer  controlled  translation  system  which  allows 
for  precise  motion  in  the  X-,  Y-  and  Z-axes.  Measurements  at  various  Y-positions  provided  radial  and  axial 
data.  In  this  paper,  ion  velocity  and  subsequently  derived  kinetic  energy  distributions  are  presented  along  the 
centerline  of  the  thruster  axis  at  positions  from  Z  =  — 6  mm  to  +40  mm  into  the  plume.  The  depth  at  which 
the  data  points  could  be  collected  inside  the  thruster  channel  was  limited  by  the  angle  inside  the  thruster 
visible  to  the  collection  optics  (60°  from  the  leading  edge  of  the  channel).  Additionally,  distributions  are 
presented  along  the  exit  plane  of  the  thruster  at  Z  =  0  mm  with  a  2  mm  spatial  resolution  in  the  Y-direction 
from  Y  =  0  mm  to  +14  mm. 


III.  Results  and  Discussion 

Two  operating  conditions  are  presented  in  this  study  of  the  ion  velocity  and  energy  distributions  of  the 
CHT-30:  a  nominal  operating  condition  and  a  condition  with  reduced  magnetic  field  strength.  The  reduced 
magnetic  field  condition  has  approximately  70  %  of  the  magnetic  field  strength,  and  provides  a  quiescent 
operating  mode  of  the  thruster.  Comparisons  are  made  between  these  two  conditions  for  positions  along 
the  centerline  and  across  the  exit  plane  of  the  thruster.  Additionally,  the  reduced  magnetic  field  case  was 
run  at  high  vacuum  chamber  background  pressure.  Comparisons  between  this  high  pressure  condition  and 
the  nominal  operating  condition  are  made  along  the  centerline  of  the  thruster.  The  velocity  and  energy 
distributions  are  each  area-normalized. 

A.  Lineshape  Analysis 

Analysis  of  the  fluorescence  measurements  made  on  the  CHT-30  provides  a  great  deal  of  information  about 
the  motion  of  the  ions  both  in  the  channel  and  in  the  plume  of  the  thruster.  A  raw  LIF  signal  is  a  convolution 
of  the  ion  velocity  distribution  function  (VDF),  and  the  transition  lineshape  (including  broadening  effects 
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Figure  4.  Coordinate  system  used  for  data  points. 
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Figure  5.  Location  of  data  points  for  axial  and  radial  LIF  measurements. 


due  to  instruments,  the  finite  linewidth  of  the  laser,  etc.).  The  transition  lineshape  is  typically  much  narrower 
than  the  VDF,17  and  therefore  the  velocity  distribution  can  often  be  represented  directly  by  the  raw  LIF 
trace.18  It  has  previously  been  shown  that  the  difference  between  a  deconvolved  VDF  and  the  raw  LIF  signal 
is  small,  less  than  20%  -  in  some  cases,  such  as  VDFs  with  low  SNR  and  broad  velocity  distributions,  it  can 
even  be  advantageous  to  forgo  deconvolution  due  to  the  additional  oscillatory  noise  it  often  produces  near 
trace  end  points.19  Therefore,  the  ion  velocity  distributions  presented  in  this  paper  are  given  by  the  raw 
fluorescence  traces. 

Figures  6  and  7  show  the  analysis  process  for  a  fluorescence  measurement  taken  at  Y  =  0  mm  and 
Z  =  20  mm.  Figure  6  shows  the  raw  absorption  and  LIF  traces  in  relative  frequency  units.  The  Doppler 
shift  relationship  is  used  between  the  stationary  absorption  reference  and  the  LIF  traces  to  transform  the 
frequency  separation  into  velocity  units,  as  shown  in  Figure  7.  In  this  figure,  the  blue  curve  represents 
the  axial  VDF,  while  the  black  curve  is  the  axial  kinetic  energy  distribution.  When  normalized  such  that 
the  total  area  under  the  curve  is  equal  to  unity,  these  traces  represent  the  probability  of  an  ion  having  a 
given  velocity  or  kinetic  energy  at  that  spatial  point.  The  kinetic  energy  is  derived  using  a  simple  1/2 rriiV2 
relationship  to  the  VDF,  where  rrii  is  the  mass  of  the  xenon  ion. 

In  Figure  7,  both  the  velocity  and  kinetic  energy  curves  are  normalized  such  that  the  peak  is  at  a  value 
of  1  (rather  than  area  normalized).  This  is  useful  in  determining  the  median  values  of  velocity  and  kinetic 
energy.  The  green  and  red  curves  represent  the  integral  of  the  velocity  and  kinetic  energy  curves,  respectively. 
The  point  where  these  integrated  lines  crosses  0.5  indicates  the  median  value  on  their  corresponding  velocity 
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Figure  6.  Example  of  LIF  and  absorption  traces  taken  at  Y  =  0  mm  and  Z  =  20  mm. 
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Figure  7.  Lineshape  analysis  for  axial  ion  velocity  and  kinetic  energy  distributions  at  Y  =  0  mm,  Z  =  20  mm. 


or  kinetic  energy  distribution.  At  Y  =  0  mm  and  Z  =  20  mm,  the  median  velocity  is  13965  m/s,  and  the 
median  kinetic  energy  is  138.47  eV.  In  most  cases  the  median  values  are  very  close  to  the  mean  values,  due 
to  the  symmetry  of  the  distributions  (i.e.  at  Y  =  0  mm  and  Z  =  20  mm,  the  mean  velocity  is  13939  m/s 
and  mean  kinetic  energy  is  138.47  eV). 

To  capture  the  range  of  velocities  or  energies  over  which  the  majority  of  ions  are  moving,  it  is  possible  to 
see  where  the  integral  curves  meet  a  certain  percentage  above  and  below  the  50%  mark.  For  one  standard 
deviation  above  and  below  the  median  value,  approximately  68%  of  the  ions  have  axial  velocities  between 
12772  m/s  and  15232  m/s  (114.62  eV  to  161.35  eV).  Similarly,  for  two  standard  deviations,  approximately 
95%  of  the  ions  have  an  axial  velocity  between  10579  m/s  and  17186  m/s  (84.32  eV  to  210.33  eV).  This 
method  for  analyzing  the  raw  fluorescence  traces  is  representative  of  the  other  data  analysis  seen  in  this 
paper. 
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B.  Centerline  Measurements 

Centerline  measurements  are  presented  for  positions  of  Z  =  —  6  mm  inside  the  thruster  channel  to 

Z  =  40  mm  into  the  plume,  as  shown  in  Figure  8. 


Figure  8. 


Location  of  data  points  for  axial  and  radial  LIF  measurements  along  the  centerline  of  the  thruster. 


1.  Ion  Velocity  Distributions 

Figure  9  shows  the  evolution  of  the  ion  velocity  distributions  for  the  CHT-30  at  a)  nominal  and  b)  reduced 
magnetic  field  conditions.  The  red,  dashed  lines  represent  the  radial  distributions,  while  the  blue  solid  lines 
represent  the  axial  distributions. 

In  the  axial  direction,  ion  velocity  distributions  are  shown  both  in  the  plume  and  at  positions  from  Z  =  - 
6  mm  to  the  exit  plane  within  the  thruster  channel.  For  the  nominal  case,  the  distributions  are  relatively 
smooth  and  narrow,  with  68%  of  the  ions  contained  within  less  than  a  4,000  m/s  range  (e.g.  2,900  m/s  to 
6,800  m/s  for  one  standard  deviation  on  each  side  of  the  median  velocity  at  Z  =  -2  mm).  For  the  reduced 
magnetic  field  condition,  the  axial  ion  velocity  distributions  are  much  broader  inside  the  thruster,  with 
less  defined  peaks  giving  ranges  of  approximately  5,000  m/s  (e.g.  470  m/s  to  5,500  m/s  for  one  standard 
deviation  from  the  median  at  z  =  -2  mm).  Comparing  the  two,  slight  acceleration  between  Z  =  -6  mm  and 
the  exit  plane  can  be  seen  in  the  nominal  case,  whereas  it  is  less  noticeable  in  the  lower  magnetic  field  case 
due  to  the  broad  distributions. 

Referring  back  to  Figure  3  which  illustrates  the  cusped  magnetic  field  profile  inside  a  CHT  -  the  magnetic 
field  lines  are  believed  to  form  equipotential  surfaces,  creating  an  electric  field  that  has  a  significant  axial 
component.2  This  causes  a  convergent  ion  flux20  to  keep  ions  away  from  the  channel  walls  and  accelerates 
the  ions  along  the  centerline.  It  therefore  makes  sense  that  the  nominal,  higher  magnetic  field  case,  would 
have  more  noticeable  acceleration  and  a  more  concentrated  velocity  distribution  inside  the  channel  due  to 
better  ion  confinement  along  the  centerline.  Overall,  however,  the  median  axial  velocities  achieved  by  the 
nominal  and  reduced  magnetic  field  cases  are  the  same  within  the  uncertainty  of  the  measurements,  at  6,500 
m/s. 

Outside  the  thruster  channel,  a  significant  portion  of  the  axial  ion  acceleration  is  seen  in  the  plume.  For 
the  nominal  case,  the  ions  achieve  axial  velocities  between  16,000  and  22,000  m/s  at  Z  =  40  mm,  with  a 
relatively  consistent  acceleration  from  the  exit  plane  through  to  the  far  field  of  the  plume.  For  the  reduced 
magnetic  field  case,  the  axial  ion  velocity  distributions  narrow  and  have  a  more  pronounced  peak  as  the  ions 
progress  further  into  the  plume. 

Figure  10  depicts  the  median  axial  ion  velocities  along  the  centerline  of  the  thruster  from  Z  =  —6  mm 
inside  the  thruster  channel  to  Z  =  40  mm  into  the  plume.  This  figure  shows  that  the  reduced  magnetic  field 
case  has  a  slightly  higher  median  axial  ion  velocity  from  Z  =  10  mm  to  Z  =  40  mm,  with  the  nominal  case 
reaching  17,200  m/s  and  the  reduced  magnetic  field  case  reaching  17,700  m/s  at  Z  =  40  mm. 

For  both  conditions,  the  radial  velocity  distribution  remains  relatively  constant  along  the  centerline  of 
the  thruster,  centered  near  0  m/s.  Previous  calculations  show  that  the  most  probable  radial  velocities  vary 
from  zero  by  no  more  than  1,200  m/s.1 
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a)  Nominal 


b)  Reduced  Magnetic  Field 
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Figure  9.  Centerline  ion  velocity  distributions  for  the  nominal  and  reduced  magnetic  field  operating  conditions. 
Red  ( - )  is  the  radial  distribution,  blue  ( — )  is  the  axial  distribution. 
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For  both  conditions,  the  radial  ion  velocities  tend  to  have  wider  distributions  at  positions  close  to  the 
exit  plane  of  the  thruster  as  compared  to  the  points  farther  into  the  plume.  This  could  be  indicative  of  the 
strongly  curved  magnetic  field  profile  seen  in  cusped-field  thrusters  that  could  produce  a  wide  distribution 
of  velocity  classes  near  the  exit  plane.  This  effect  diminishes  in  strength  farther  downstream,  due  to  the 
reduced  influence  of  the  radial  component  of  electric  field  farther  from  the  exit  plane. 

Alternatively,  the  field  strength  may  be  strong  enough  to  produce  Zeeman  broadening.  However,  by 
comparing  the  two  operating  conditions,  it  is  noticeable  that  the  reduced  magnetic  field  case  tends  to  have 
wider  radial  ion  velocity  distributions,  but  sharper  peaks  in  the  velocity  than  the  nominal  case  -  indicating 
that  while  Zeeman  broadening  may  be  present,  it  is  not  the  dominant  broadening  mechanism.  Near  the  exit 
plane  of  the  thruster,  the  wide  radial  distributions  in  the  lower  magnetic  field  case  could  be  a  function  of 
the  ions  with  slower  axial  velocities  being  pulled  more  effectively  in  the  radial  direction  by  the  strong  radial 
component  of  electric  field.  Similar  to  the  axial  velocity  distributions,  a  stronger  magnetic  field  seems  to 
produce  less  variation  in  radial  velocities. 


Figure  10.  Median  axial  velocities  along  centerline  of  thruster  plume  for  nominal  and  reduced  magnetic  field 
conditions. 


2.  Ion  Kinetic  Energy  Distributions 

Figure  11  shows  the  ion  kinetic  energy  distributions  at  the  centerline  of  the  thruster  channel  for  the  nominal 
(red)  and  reduced  magnetic  field  (blue)  conditions.  Previous  studies  on  the  CHT  showed  that  the  majority 
of  the  xenon  propellant  is  ionized  at  the  border  of  the  annular  and  cylindrical  region  inside  the  CHT-30. 
The  ion  acceleration  then  occurs  in  the  cylindrical  region  of  the  channel,  along  which  the  greatest  potential 
drop  occurs.20  However,  comparing  the  centerline  kinetic  energies  achieved  by  the  exit  plane  of  the  thruster 
channel  to  those  achieved  outside  shows  that  a  large  portion  of  the  axial  acceleration  actually  occurs  in  the 
thruster  plume,  at  positions  further  than  Z  =  5  mm. 

Throughout  the  plume,  the  width  of  the  kinetic  energy  distributions  remains  similar  for  both  operating 
conditions.  However,  at  Z  =  40  mm,  the  distribution  for  the  nominal  condition  becomes  much  wider  than 
that  of  the  reduced  magnetic  field  case  (i.e.  the  full- width  at  half  maximum  (FWHM)  is  25  eV  for  the  reduced 
magnetic  field  case,  and  50  eV  for  the  nominal  case).  The  wings  of  the  energy  distributions,  however,  when 
measured  by  one  standard  deviation  on  either  side  of  the  median  kinetic  energy  value  give  a  range  of  177  eV 
to  228  eV  for  the  nominal  case,  and  191  eV  to  270  eV  for  the  reduced  magnetic  field  case  -  a  wider  range 
of  higher  energy  ions  for  the  off-nominal  case.  The  reason  the  range  over  which  the  reduced  magnetic  field 
kinetic  energies  is  wider  in  terms  of  percentage  of  ions  is  that  while  the  peak  in  the  kinetic  energy  distribution 
is  sharper  for  this  condition,  the  wings  are  much  shallower,  favoring  the  higher  energy  side  of  the  distribution. 
The  reduced  magnetic  field  condition  appears  to  accelerate  the  ions  more  effectively  at  further  distances  in 
the  plume. 

Figure  12  shows  the  median  axial  kinetic  energies  along  the  centerline  of  the  plume.  The  median  value 
for  both  the  nominal  and  reduced  magnetic  field  kinetic  energy  is  approximately  25  eV  at  the  exit  plane  of 
the  thruster,  with  an  additional  175  eV  to  188  eV  gain  in  the  plume.  At  a  position  40  mm  into  the  plume, 
this  constitutes  an  energy  of  200  eV  for  the  nominal  case  and  213  eV  for  the  reduced  magnetic  field  case 
for  an  applied  potential  of  300  V.  This  confirms  that  the  acceleration  is  not,  in  fact,  primarily  inside  the 
thruster  channel. 
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Figure  11.  Centerline  ion  kinetic  energy  distributions  for  the  nominal  operating  condition  in  red  and  reduced 
magnetic  field  operating  condition  in  blue. 
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Figure  12.  Median  axial  kinetic  energies  along  centerline  of  plume  for  nominal  and  reduced  magnetic  field 
conditions. 
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Figure  13.  Centerline  axial  electric  fields  derived  from  median  kinetic  energies  for  nominal  and  reduced 
magnetic  field  conditions. 

Figure  13  shows  the  axial  electric  field  along  the  centerline  of  the  thruster,  derived  from  the  kinetic  energy 
values.  Comparing  the  electric  fields  makes  the  differences  in  acceleration  seen  in  the  velocity  and  kinetic 
energy  distributions  more  clear.  The  reduced  magnetic  field  case  has  a  noticeably  higher  axial  electric  field 
around  Z  =  10  mm  to  20  mm,  which  explains  how  the  ions  gain  more  energy  in  this  region  as  compared  to 
the  nominal  condition. 

C.  Centerline  Measurements  at  High  Pressure 

The  reduced  magnetic  field  condition  was  also  run  with  a  higher  chamber  background  pressure  of  3.2  x 
10-5  torr,  for  several  positions  throughout  the  plume.  Figure  14  shows  the  velocity  distributions  in  the 
radial  and  axial  directions  along  the  centerline  of  the  thruster.  As  seen  in  the  previous  conditions,  the  radial 
velocity  distributions  remain  relatively  constant,  centered  around  zero  throughout  the  plume.  For  the  axial 
direction,  however,  there  is  a  large  increase  in  velocity  between  Z  —  0  mm  and  5  mm,  not  seen  in  the  lower 
pressure  conditions.  This  is  easily  seen  in  Figure  15  which  shows  the  median  velocities  along  the  centerline 
of  the  thruster. 

Figure  16  shows  the  axial  kinetic  energy  distributions  for  this  higher  pressure  case,  with  the  nominal 
condition  also  shown  for  comparison.  It  is  evident  that  between  Z  =  0  mm  and  5  mm,  the  gain  in  kinetic 
energy  by  far  surpasses  that  seen  in  the  nominal  condition.  But  by  Z  =  40  mm,  the  axial  ion  kinetic  energies 
become  nearly  the  same. 

This  is  summarized  in  Figure  17,  which  shows  the  median  kinetic  energies  along  the  centerline  of  the 
thruster.  For  both  the  nominal  and  high  pressure  reduced  magnetic  field  cases,  the  kinetic  energy  is  around 
25  eV  at  the  exit  plane  of  the  thruster.  The  high  pressure  case  gains  nearly  100  eV  of  energy  in  the  first  10 
mm  of  the  plume  whereas  the  nominal  case  gains  approximately  20  eV.  By  Z  =  40  mm,  the  ion  energies  are 
once  again  very  comparable  with  median  values  near  200  eV. 


12  of  20 


American  Institute  of  Aeronautics  and  Astronautics 
Approved  for  public  release;  distribution  is  unlimited. 


Reduced  Magnetic  Field  at  High  Pressure 
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Figure  14.  Effect  of  high  pressure  on  reduced  magnetic  field  condition  velocity  distributions.  Red  ( - )  is 

the  radial  distribution,  blue  ( — )  is  the  axial  distribution. 
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Figure  15.  Median  axial  velocities  along  centerline  of  thruster  plume  for  nominal  and  high  pressure  reduced 
magnetic  field  conditions. 


13  of  20 


American  Institute  of  Aeronautics  and  Astronautics 
Approved  for  public  release;  distribution  is  unlimited. 


Kinetic  Energy  [eV] 

Figure  16.  Effect  of  high  pressure  on  reduced  magnetic  field  condition  kinetic  energy  distributiosn.  Blue  is 

the  reduced  magnetic  field  at  high  pressure  condition,  while  red  ( - )  is  the  nominal  condition,  shown  for 

reference. 


Figure  17.  Median  axial  kinetic  energies  along  centerline  of  thruster  for  nominal  and  high  pressure  reduced 
magnetic  field  conditions. 
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Figure  18.  Centerline  axial  electric  fields  derived  from  median  kinetic  energies  for  nominal  and  high  pressure 
reduced  magnetic  field  conditions. 

For  both  cases,  the  wings  of  the  velocity  distributions  stretch  between  approximately  175  eV  and  225  eV 
for  68%  of  the  ions,  showing  that  a  good  portion  of  the  applied  300  V  anode  potential  is  used  to  accelerate 
the  ions  in  the  axial  direction,  but  still  less  than  in  the  reduced  magnetic  field  case  at  lower  pressure. 

Figure  18  compares  the  electric  field  of  the  reduced  magnetic  field  at  high  pressure  condition  to  the 
nominal  condition.  This  figure  shows  very  distinctly  that  raising  the  background  pressure  in  a  vacuum 
chamber  moves  the  acceleration  region  upstream  towards  the  exit  plane  of  the  thruster.  This  effect  has  been 
seen  in  previous  studies.21 

D.  Exit-plane  Measurements 

Exit  plane  measurements  are  presented  for  positions  of  Y  =  0  mm  to  14  mm  along  the  face  of  the  thruster 
channel  at  Z  =  0  mm,  as  shown  in  Figure  19. 
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Figure  19.  Location  of  data  points  for  axial  and  radial  LIF  measurements  along  the  centerline. 


1.  Ion  Velocity  Distributions 

Figure  20  shows  the  progression  of  ion  velocity  distributions  across  the  exit  plane  of  the  thruster  channel. 
Both  the  nominal  and  reduced  magnetic  field  conditions  behave  similarly  across  the  exit  plane.  The  axial 
velocities  slowly  decrease  across  the  width  of  the  channel  towards  Y  =  14  mm.  The  axial  distributions 
broaden  at  positions  farther  from  the  centerline,  indicating  a  wider  range  of  ion  velocities  -  in  particular,  a 
larger  percentage  of  ions  with  slow  axial  velocities. 
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a)  Nominal 


b)  Reduced  Magnetic  Field 
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Figure  20.  Exit  plane  ion  velocity  distributions  for  the  nominal  and  reduced  magnetic  field  operating  conditions. 
Red  ( - )  is  the  radial  distribution,  blue  ( — )  is  the  axial  distribution. 
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For  example,  the  nominal  condition  at  Y  =  0  mm  has  a  median  axial  velocity  of  5960  m/s,  with  wings 
spanning  3160  m/s  to  7870  m/s  for  68%  of  the  ions.  At  Y  =  14  mm,  the  median  axial  velocity  is  slower  at 
4630  m/s,  with  a  wider  span  of  1850  m/s  to  8640  m/s.  The  reduced  magnetic  field  case  loses  slightly  less 
axial  velocity,  starting  at  6600  m/s  at  7  =  0  mm  and  decreasing  to  5860  m/s  at  the  channel  wall.  For  both 
cases,  this  could  be  reminiscent  of  the  convergent  ion  flux  inside  the  thruster  channel  previously  described 
by  Smirnov,  et  ah, 20  where  the  more  energetic  ions  are  kept  away  from  channel  walls. 

Also  of  note  in  the  exit  plane  ion  velocity  evolutions  is  the  increase  in  radial  velocity  towards  the  channel 
wall.  Starting  at  the  centerline,  the  radial  component  of  ion  velocity  is  nearly  zero  for  both  operating 
conditions.  But  at  a  position  of  Y  =  14  mm,  the  radial  velocities  become  comparable  in  magnitude  to 
the  axial  velocities.  Figure  21  gives  the  median  values  -  4360  m/s  for  the  nominal  condition  (compared 
to  4630  m/s  in  the  axial  direction)  and  4870  m/s  for  the  reduced  magnetic  field  condition  (compared  to 
5860  m/s  in  the  axial  direction).  This  high  radial  component  of  ion  velocity  at  the  exit  plane  is  indicative  of 
the  large  plume  divergence  half  angle  that  is  typical  of  cylindrical  Hall  thrusters,  but  is  very  different  than 
traditional  annular  Hall  thrusters  where  the  axial  component  of  velocity  is  much  higher  than  the  radial.12 


Y-Position  [mm] 

Figure  21.  Median  radial  velocities  along  exit  plane  of  thruster  for  nominal  and  reduced  magnetic  field 
conditions. 


2.  Ion  Kinetic  Energy  Distributions 

Figure  22  shows  the  radial  kinetic  energy  distributions  across  the  exit  plane  of  the  thruster.  For  both 
operating  conditions,  between  positions  of  Y  =  0  mm  and  Y  =  12  mm,  the  peaks  in  the  distributions 
are  centered  very  close  to  0  eV.  Towards  the  channel  wall,  however,  the  tails  of  the  distributions  begin  to 
stretch  towards  higher  radial  energies.  The  median  radial  kinetic  energies  are  summarized  in  Figure  23.  At 
Y  —  14  mm,  the  median  radial  kinetic  energy  is  12  eV  for  the  nominal  case  and  14  eV  for  the  reduced 
magnetic  field  case.  Adding  these  values  to  the  nominal  axial  kinetic  energies,  15  eV  for  the  nominal  case 
and  23  eV  for  the  reduced  magnetic  field  case  at  Y  =  14  mm,  shows  that  there  is  actually  a  gain  in  energy 
towards  the  channel  wall  at  the  exit  plane.  Both  conditions  start  at  approximately  25  eV  at  Y  =  0  mm. 
The  nominal  condition  increases  to  26  eV  at  the  channel  wall,  and  the  reduced  magnetic  field  case  increases 
to  37  eV.  Overall,  while  there  appears  to  be  a  convergent  ion  flux  inside  the  thruster  channel,  the  kinetic 
energy  distributions  at  the  exit  plane  indicate  that  the  ions  are  very  divergent  just  outside  the  channel. 

Figure  24  shows  the  radial  electric  field  profile  across  the  exit  plane  of  the  thruster.  Previous  studies 
have  shown  the  half  angle  of  the  plume  to  be  as  large  as  60°  to  80°,  compared  to  traditional  annular  Hall 
thrusters  which  have  plume  divergence  half  angles  of  approximately  45°  to  50°. 4  This  increased  divergence 
angle  could  be  due  to  the  radial  component  of  the  electric  field  in  this  region  being  stronger  than  that  of  a 
traditional  annular  Hall  thruster.2  This  can  be  explained  by  the  evolution  of  the  magnetic  field  profile  as  it 
gets  farther  from  the  centerline  of  the  channel.  The  resulting  equipotential  surfaces  provide  an  electric  field 
profile  with  a  significant  outward  pointing  radial  component.  From  the  exit  plane  outwards  into  the  plume, 
all  radial  velocity  components  are  therefore  pointing  away  from  the  centerline.  This  plays  a  large  role  in 
the  divergence  of  the  plume,  thereby  reducing  or  eliminating  the  centerline  ion  flux  seen  inside  the  channel, 
resulting  in  a  reduced  mixing  zone  compared  to  annular  Hall  thrusters. 
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Y  =  6  mm 


Y  =  10  mm 


Y  =  12  mm 


Kinetic  Energy  [eV] 


Figure  22.  Median  radial  kinetic  energy  distributions  along  exit  plane  of  thruster  for  nominal  and  reduced 
magnetic  field  conditions. 
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Y-Position  [mm] 

Figure  23.  Median  radial  kinetic  energies  along  exit  plane  of  thruster  for  nominal  and  reduced  magnetic  field 
conditions. 


Y-Position  [mm] 

Figure  24.  Exit  plane  radial  electric  fields  derived  from  kinetic  energies  for  nominal  and  reduced  magnetic 
field  conditions. 


IV.  Conclusions 

Laser  induced  fluorescence  was  used  to  map  the  evolution  of  the  ion  velocity  and  kinetic  energy  distri¬ 
butions  of  the  CHT-30  both  inside  the  thruster  channel  and  in  the  plume.  A  nominal  operating  condition 
was  examined  and  compared  to  an  reduced  magnetic  field  case  and  the  reduced  magnetic  field  case  at  an 
elevated  background  pressure. 

Inside  the  thruster  channel,  reduction  of  the  magnetic  field  strength  created  wide  velocity  distributions 
with  a  larger  portion  of  ions  with  slower  axial  velocities  near  the  exit  plane  of  the  thruster.  The  nominal 
condition  appeared  to  have  more  axial  ion  acceleration  inside  the  channel,  while  the  reduced  magnetic  field 
case  had  a  greater  acceleration  in  the  plume,  especially  after  Z  m  10  mm.  By  40  mm  into  the  plume,  the 
median  axial  ion  kinetic  energies  were  200  eV  for  the  nominal  condition  and  213  eV  for  the  reduced  magnetic 
field  condition,  with  wings  of  the  energy  distributions  approaching  270  eV,  showing  that  a  large  portion  of 
the  applied  300  V  anode  potential  went  mostly  into  axially  accelerating  the  ions. 

Raising  the  background  pressure  had  the  effect  of  pushing  the  axial  acceleration  upstream  towards  the  exit 
plane.  This  condition  had  the  majority  of  its  acceleration  around  Z  =  5  mm.  By  Z  =  40  mm,  the  velocities 
in  the  nominal  condition  caught  up  with  the  high  pressure  condition,  resulting  in  similar  distributions.  The 
median  axial  kinetic  energy  for  the  high  pressure  case  was  approximately  200  eV  at  Z  =  40  mm,  with  68% 
of  the  ions  having  energies  between  175  eV  and  225  eV. 

Looking  at  the  distributions  across  the  exit  plane  of  the  thruster,  for  both  the  nominal  and  reduced 
magnetic  field  operating  conditions  the  radial  component  of  ion  velocity  started  near  zero,  then  increased  to 
be  comparable  in  magnitude  to  the  axial  velocities  near  the  channel  wall.  This  agreed  well  with  the  large 
plume  divergence  angles  seen  in  these  type  of  thrusters,  and  indicated  that  the  convergent  ion  flux  seen 
inside  the  thruster  channel  disappears  into  a  dispersed  plume. 

Along  with  an  increased  axial  acceleration  along  the  centerline  of  the  plume,  the  increased  radial  velocity 
in  the  reduced  magnetic  field  case  implied  that  reducing  the  current  to  the  electromagnets  not  only  changed 
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the  magnetic  field  strength,  but  also  the  topology.  Future  work  includes  comparison  of  these  CHT-30  results 
to  LIF  measurements  on  other  types  of  cusped  field  thrusters,  such  as  the  MIT  and  Stanford  Diverging 
Cusped  Field  Thrusters  (DCFTs). 
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